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ABSTRACT 

I re-examine the relation between virial mass and concentration for groups and clusters of 
galaxies as measured in a number of recent works. As previously noted by several authors, 
low-mass clusters and groups of galaxies display systematically larger concentrations than 
simple prescriptions based on pure n-body simulations would predict. This implies an ob- 
served concentration-mass relation with a substantially larger slope/normalization than ex- 
pected from theoretical investigations. Additionally, this conclusion seems to be quite inde- 
pendent on selection effects, holding for both lensing based and X-ray based cluster samples. 
In order to shed new light on this issue I employ a simple spherical halo model containing, 
in addition to dark matter, also stars and hot diffuse gas in proportions and with distributions 
in agreement with the most recent observations. Moreover, I include the contraction effect 
experienced by dark matter due to the cooling of baryons in the very central part of the struc- 
ture itself. The resulting modified concentration-mass relation is steeper than the theoretical 
input one, because star formation is fractionally more efficient in low-mass objects. How- 
ever, the effect is non-vanishing at all masses, thus resulting also in a larger normalization. 
Overall the new relation provides a better representation of the observed one for almost all 
catalogs considered in this work, although the specific details depend quite significantly on 
the baryon fraction prescription adopted. Specifically, the observed concentration-mass rela- 
tion seems to favor a scenario where the stellar mass fraction in large clusters of galaxies is 
substantially lower than several works have found. Anyhow, the same effect could also be 
produced by a redistribution of baryons within the structure. Moreover, the concentration of a 
number of high-mass objects seems to be significantly lower even than the predictions based 
on pure n-body simulations, and they are hence unaccounted for in the modified scenario that 
is proposed here. Finally I use this simple model to show how the estimated concentration of 
cosmic structures is expected to be overestimated as a function of the radial range covered by 
the analysis. 
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1 INTRODUCTION 

One of the out standing successes of the standard ACD M cosmo- 
logical model dRubin et al.l |2009| ; iKomatsu et al.l 1201 lb is its ac- 
count for the formation of structures in the Universe. Accordingly, 
tiny density fluctuations were laid down in the dark matter distri- 
butio n during the early stages of the Universe (Hawking & Moss! 
Il983l) . These seed perturbations grew up in time due to gravita- 
tional instability, and eventually detached away from the overall 
expansion and collapsed, giving rise to the web of clumps, fila- 
ments, and vo ids that compose th e cur rently observed large scale 
structure (e.g JPercival et alj20oTll2007h . 

The redshift evolution of dark matter structures is determined 
only by gravity, and it is hence quite easy to model with the help 
of numerical n-body simulations. As a matter of fact, simulations 
provide a plethora of predictions that have been tested against 
observations, up to now with an overall good degree of success. 
One of the principal predictions is that dark matter bound struc- 



tures (or halos) should have a cuspy d ensity profile, on average 
well r epresented by the two-parameter Navarro, Frenk, & White 
( 1996) (NFW henceforth, see also iDubinski & Carlberd 1 199lF 
iNavarro. Frenk. & White! 1995Lfl997h function 



p(x) : 



Ps 



(1) 



x(l +x) 2 ' 

where p s is a scale density, while x = r/r s is the radial distance 
from the center of the halo in units of a scale radius r s . The two pa- 
rameters of the NFW profile can be easily related to more directly 
measurable quantities, namely the mass and the concentration of 
the halo. The mass of a dark matter halo is conventionally defined 
as the mass contained within a sphere whose average density equals 
some factor A times the critical density of the Universe at the red- 
shift of interest, 



m A = -nR\ [Ap c (z)] 



(2) 



The specific value of A varies according to the author. The most 
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common approach is to set A = 200, independently on redshift and 
cosmology. It has been shown that the R 2 m thus defined is a good 
approximation of the radius within witch the s tructure is in dynami- 
cal equilibrium dEke. Navarro. & Frenkl 19980 . Another popular ap- 
proach is to have A = A v (z), where A v (z) is the virial overdensity 
according to the sphe rical collapse model, wh ich depends on red- 
shift and cosmology (Bry an & Normanll 19981) . Also, a number of 
authors prefer to replace the critical density of the Universe with 
the mean matter density. In this work I adopted the former choice, 
with A = 200. The concentration of a halo is a measure of its com- 
pactness, and is defined for a NFW profile as ca = /?a/jV m what 
follows I will drop the suffix A in all quantities, leaving implicit 
that A = 200, and I shall refer torn = mioo and R = Z?20o as to virial 
quantities. 

Another cornerstone prediction of the standard cosmology is 
the hierarchical clustering, meaning the fact that high-mass objects 
are formed as aggregates of lower-mass ones, that were formed at 
earlier times in the same manner. A related feature is that the con- 
centration and the mass of a dark matter structure are not unre- 
lated, rather the concentration is a decreasing function of mass. 
This result is commonly interpreted as due to the fact that ha- 
los retain information on the average matter density of the Uni- 
verse at the moment of their formation dNavarro. Frenk. & White! 
1996). Since low-mass structures form earlier than large-mass ones, 
they have more time to relax and compactify, hence displaying 
a larger concentration. The concentration-mass relation predicted 
by n-body simulations is well represe nted by a shallow power law 
c(m,z = 0) oc m-", wi th a ~ 0.1 dDolag et al.ll2004 IGao et al.1 
l2008t IZhaoetalJl2009T) . Moreover, the distribution around this 
mean relation is basically log-normal, with a logarithmic mean de- 
viation ranging from ~ 0. 15 for th e most rela xed structures, up to 
~ 0.30 for the most disturbed ones ( Jing 2000). 

The redshift dependence of the concentration for a given mass 
is much less established. While the original NFW prescription dis- 
plays v ery little redshift evolution, simple arguments based o n later 
studies dBullock et al.l200ll ; lEke. Navarro. & Steinmetj200lh sug- 
gest that the concentration of a dark matter halo with fixed mass 
should evolve quite strongly with redshift, basically oc (1 + zT 1 ■ 
More recent studies however dNeto et alj|2007l; iDuffv et alj|2008l ; 
IGao et al.ll2008l ; IZhao et alJl2009l ; |Prada et al.ll201lh imply that the 
original NFW prediction was not so wrong after all. The concen- 
tration of high-mass objects shows basically no evolution with red- 
shift, while the overall normalization evolves much less strongly 
than oc (1 + z)~\ as previously thought. 

Substantial complication is added to this picture by the role of 
baryonic matter, that is gas and stars. Although these luminous mat- 
ter components are almost always subdominant with respect to dark 
matter, and hence are not expected to drive the process of structure 
formation, they are well known to have an impact on the details 
of matter allotment dJingetalJl2006l : lFedeli. Dola gT& Moscardinil 
|20T 1). While on very large scales (typically above a few tens of 
Mpc) it is safe to assume that the gas distribution follows the dark 
matter one, on small scales (say below 0. 1 Mpc) the gas becomes 
dense enough to cool down radiatively and condense into stars. 
Massive stars themselves subsequently provide energy and metal 
injection into the gas when they go off as supernovae, with Ac- 
tive Galactic Nuclei (AGN) activity constituting another important 
source of energetic feedback. All these non-gravitational processes 
are very difficult to model, and the majority of them takes place 
at scales that are too small to be resolved by current cosmological 
simulations, thus requiring the implementation of sub-grid approx- 
imations. 



As of today, no consensus has been reached as for the impact 
of baryons on the small scale matter distribution. Yet, observations 
are now achieving a precision level where such an impact cannot 
be neglected anymore. Examples in this direction are given by fu- 
ture wide field weak lensing surveys (such as Euclid and WFIRST, 
lLaureiis et al.l |20 fib that will be capable of measuring the non- 
linear matter power spe ctrum at the percent l evel, or the upcoming 
HST CLASH program (Postma n et alj[2oTTb . that will determine 
the internal structure of massive galaxy clusters with a relative ac- 
curacy of ~ 5%. In both cases the impact of baryons is expected to 
be substantially larger than the forecasted measurement precisions. 

In this paper I used a spherical semi-analytic model in or- 
der to investigate the impact of baryonic physics, specifically gas 
cooling and star formation, on the concentration-mass relation of 
galaxy systems. I built upon several observational results, showing 
that low-mass clusters and groups of galaxies are systematically 
over-concentrated with respect to theoretical predictions. When- 
ever possible, I adopted observationally motivated prescriptions as 
ingredients of the model, in order to partially bypass the uncer- 
tainty typical of numerical simulations stemming from the amount 
of cooling/feedback that is requested to reproduce a wide array 
of obs ervable results. In a recent alternative approach, Duf fy et all 
d2010h used a suite of cosmological simulations where different 
kinds of baryonic physics were implemented in order to study their 
impact on the structure of groups and clusters of galaxies. Dur- 
ing the course of the paper I shall compare my results to those of 
Duff y et all d20ld) where relevant. 

The rest of this work is organized as follows. In Section [2] I 
show the concentration-mass relation recently observed in differ- 
ent cluster catalogs by several authors. In order to allow a self- 
consistent comparison I accurately convert all masses and concen- 
trations to the convention adopted in this work. In Section [3] I de- 
scribe in detail the various ingredients entering the simple semi- 
analytic model that I employed in order to obtain a modified the- 
oretical concentration-mass relation. In Sectional summarize the 
results, and in Sectional discuss them at length. Section|6]is de- 
voted to my conclusions. In what follows I adopted the standard flat 
ACDM cosmological model, with parameter values suggested by 
the latest WMAP data analysis dKomatsu et al.l20Tlh . This implies 
a matter density parameter O m = 0.272, a cosmological constant 
density parameter f2 A ,o = 1 - fl m ,o, a Hubble constant Hq = h 100 
km s -1 Mpc -1 , with h = 0.704, and a normalization of the matter 
power spectrum given by cr 8 = 0.809. 



2 OBSERVED CONCENTRATIONS 

In this Section I describe the concentration-mass relation resulting 
from several different observational studies of groups and clusters 
of galaxies. As mentioned in Section [T] different authors in gen- 
eral adopt different conventions for defining the size (and hence the 
mass and concentration) of a structure, and in some occasion even 
slightly different cosmological parameter values. In order to make 
self-consistent comparisons it is hence necessary to convert all val- 
ues of mass and concentration to the cosmology and convention 
adopted in this work. As a reminder, I define the size of a structure 
as the radius of the sphere whose average density equals to Ap c (z), 
where A = 200, and p c (z) = 3H 2 (z)/8nG is the critical density of 
the Universe at the redshift of the selected cluster. Please note that 
the critical density of the Universe and the average matter density 
are related via the matter density parameter f2 m (z), thus choosing 
the latter instead of the former as a reference density ultimately 
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translates in a different value of A as a function of redshift. Analo- 
gously, different cosmologies produce a different p c (z), which can 
also be regarded as an altered overdensity value as a function of 
redshift. 

The procedure for converting masses and concentrations mea- 
sured with respect to a certain overdensity A[ into those measured 
with respect to another overdensity A 2 by assuming a NFW profile 
is relatively stra i ghtfor ward, and it is condensed in the Appendix of 
Hu & Kravtsov (2003). Summarizing very briefly, the correspond- 
ing concentrations can be linked by numerically solving the alge- 
braic equation 



F(ci) A! F(c 2 ) 



0. 



(3) 



(for which lHu & Kravtsovi 2003 give an approximate analytical so- 
lution) while the corresponding masses are linked by 

\3 A 



The function F{x) reads 
F(x) = ln(l + x) - 



(4) 



(5) 



1 + X ' 

and can be used to relate the mass with the concentration of a NFW 
dark matter halo according to M = 4np s rlF(c). It is easy to see that 
if A] > A 2 , then C\ < c 2 , and also Mj < M 2 , because the radius of 
the structure must encompass a larger mean density. 

In numerous of the observational works described below the 
authors refer to a modified concentration c (m, z) = (1 + z) c(m, z), 
rather than to the original concentration. The reason for this is the 
fact mentioned in Section[T]that according to some old theoretical 
concentration-mass prescriptions the concentration at a fixed mass 
is expected to scale with redshift as oc ( 1 +z)~ 1 . Thus, mutiplying the 
concentrations by (1 + z) would arguably remove this redshift de- 
pendence, and hence render measurements for objects at different 
redshifts more directly comparable. However, several more recent 
numerical studies have shown that the redshift dependence of the 
concentration is much weaker than that, and likely even vanishing 
for high-mass clusters. This has been verified observationally by 
ISchmidt & Alien! d2007h (see below for more details) who found 
a redshift dependence of the concentration-mass normalization of 
oc (1 + z )-°-7i±o,52_ j t j s eaS y to understand that assuming a redshift 
dependence of the concentration significantly stronger than the ac- 
tual one and then correcting for it would bias high the concentration 
of the highest redshift objects, that are presumably the most mas- 
sive ones due to selection effects (especially in X-ray studies). This 
would have the effect of making the c - M relation so mewhat flat- 
ter tha n it actually is. Given this, and in agreement with Og uri et all 
I chose not to include any kind of redshift correction, thus 
assuming that all the measured concentrations are as they would be 
at z = 0. The weak redshift dependence of the concentration and 
the limited redshift range of the majority of the studies described 
below suggest that this is a good approximation, however I shall 
discuss more on it in Section[5] 

2.1 The data 

I considered six of the most recent compilations of galaxy groups 
and clusters for the purpose of measuring the concentration-mass 
relation. In the following I describe each of them in detail, and 
summarize their findings. In Figure[T]l show the position (with rel- 
ative errors) of objects in each catalog in the mass-concentration 



plane, as well a s the expectation from the theoretical study of 
IGao et al. I J2008h and the best fitting power-law for the observed 
clusters. Specifically, I fit the concentration-mass relation with the 
two-parameter function 



(in 
— 
m 



(6) 



where mo = 5 x 10 14 h~' M , and the concentration for a fixed 
mass is assumed to be redshift-independent. In Figure Q] the val- 
ues of mass and concentration for each cluster have been carefully 
converted to the overdensity convention and fiducial cosmological 
model adopted in the present paper, as discussed above. 

The observational data are fitted with the power law of Eq. 
((6), and the best fitting parameters are found by minimizing the 
X 2 ( c o, a ) function defined as 



X 2 (c ,a) = 



[Log a - Log c{mu z = 0\ cq, a)] 2 



(7) 



where n is the number of galaxy systems in the catalog, <x, is the 
logarithmic error on the measurement of the z-th concentration C;, 
and & is the expected intrinsic scatter in the concentration at a given 
mass, that we set to & - 0.15 (see Section[T). One should bear in 
mind that the best fitting normalization and slope depend on the 
chosen pivot mass, and that for an arbitrarily chosen mo the two 
shall have some degree of correlation. Nonetheless I am not inter- 
ested here in finding the best possible fit to the data, rather only to 
demonstrate the degree of disagreement of the latter with theoreti- 
cal models. 



2.7.7 \Buoteet ali\200^) 

In this work the authors used a sample of 39 galactic systems 
spanning a wide range of masses, from isolated ellipticals (m ~ 
5 x 10 12 h' l M & ) to massive galaxy clusters (m > 10 15 /r'M G ). 
The objects were chosen to have the best quality Chandra and 
XMM-Newton observations to date. Moreover, in order to ensure 
hydrostatic equilibrium of the hot gas to be a good approximation, 
only systems displaying very regular X-ray maps, devoid of strong 
asymmetries, were selected. All objects in the sample are relatively 
local, having redshifts z < 0.2 and with only a handful having 
z > OA, 

IBuote et al. (2007) used the aforementioned X-ray observa- 
tions to infer the spherically averaged mass profile of the objects 
in the sample, and hence determined the virial mass and the con- 
centration for each one. It is worth stressing that in their anal- 
ysis the authors corrected the concentrations by the extra factor 
(1 + z). As noted above, this possibly introduces a bias, making the 
concentration-mass relation flatter than it actually is. Despite this, 
they still found a slope of a = 0.17 , steeper than the th eoretical ex- 
pectation of a ~ 0.1. Additionally. IBuote et aljj2007t) came across 
the fact that low-mass objects tend to lie systematically above the 
best-fitting relation, and that by including only massive systems in 
their statistical analysis the slope of the concentration-mass rela- 
tion decreases substantially. These findings suggest that low-mass 
structure are responsible for the steepness of the relation. 

The authors also discuss the normalization of the 
concentration-mass relation, explaining that it is consistent 
with the expectations for a ACDM cosmology having a high matter 
power spectrum normalization <x 8 = 0.9, although as explained 
previously low mass systems are systematically above that. This, 
and the fact that recent observations favor a erg value actually 
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closer to 0.8 rather than 0.9 dKomatsu et alj2oTlh . would advocate 
again for some mechanism to increase halo concentrations with 
respect to the predictions of pure n-body simulations. The top 
left panel of Figure Q] summarizes the situation for this catalog. 
Without correcting for any redshift dependence, I obtain a best fit 
slope of a = 0.20 ± 0.04, margina lly steeper than the one obtained 
originally by Buote et alj (2007). The theoretical prediction by 
iGao et al U2OO8F is significantly flatter than required by the data, 
with only 4 out of 39 observed structures lying below it. 



2.1.2 \Comerford & Nataraiaii feoO'A) 

In this work the authors presented a compilation of virial mass 
and concentration measurements derived in different works and us- 
ing different techniques. Specifically, the mass distributions of 100 
galaxy clusters and groups were obtained using X-ray observations, 
gravitational lensing (weak and/or strong), kinematic analysis, i.e., 
line-of-sight velocity dispersions, and the caustic method. The fact 
that different techniques are mixed together is a significant variant 
with respect to other works, where commonly all cluster parame- 
ters are measured adopting the same technique. This makes it more 
difficult to interpret the results, as different methods are (or may be) 
affected by different biases. The cluster sample spans a mass range 
included between in ~ 3 x 10 13 hr l M and m ~ 2x 10' 5 /r'M , 
and comprises very local systems as well as a few distant ones (up 
to z ~ 0.9). 

Similarly to iBuote et ID d2007h . the authors multiplied the 
cluster concentrations by the extra factor (1 + z), thus possibly 
biasing the slope of the concentration-mass relation low. Accord- 
ing to t heir analysis , this slope measures a = 0.15, consistent 
with the lBuote et al.1 {2007) result and marginally steeper than the 
theoretical expectation. Also, the normalization of the relation is 
substantially higher than, and basically in disagreement with, the 
one predicted for a cosmology with <x 8 ~ 0.8. In the top right 
panel of Figure [T] I report the concentration-mass r elation for the 
catalog collected by Comerford & Nataraja^ d2007l) . Similarly to 
what they did, I also removed all objects that have no mass mea- 
surements or no reported errors for either the mass or the con- 
centration. When more than one measurement is available for the 
same cluster, I use the average values, with errors equal to the 
largest errors reported. The observed scatter is very large, result- 
ing in a slope of a = 0.12 ± . 06, co nsistent with both the the- 
or etical prediction of Gao et al. I J2008h and the original analysis 
of IComerford & NataraiarJ J2007I) . The normalization on the other 
hand is substantially larger than expected for a WMAP-based cos- 
mological model. 



2.1.3 \Ettoriet alU201(\) 

These authors used X-ray surface brightness and temperature pro- 
files measured by XMM-Newton for a sample of 44 luminous 
galaxy clusters at 0. 1 < z < 0.3 to infer their mass distributions, 
and h ence the re l ation between the virial mass and the concentra- 
tion. lEttori et aL U20ld) excluded from their sample objects show- 
ing signs of recent and strong interactions, in order to ensure the 
hydrostatic equilibrium to be a good approximation for the hot gas. 

They used two different methods in order to reconstruct the 
mass profiles of clusters in their samples, which give consistent re- 
sults between each other. I shall refer to the results of their first 
method, as it has been extensively tested against numerical simula- 
tions. In this work as well the authors corrected each concentration 



value by the factor (1 + z), and found that both the slope and the 
normalization of the relation agree with expectations from numer- 
ical simulations. However, they also claim that the range in mass 
they explore, m > 10 14 fr 1 M is too small to draw any definitive 
conclusion. 

In the middle left panel of Figure[T]l repo rt the concentration - 
mass relation as measured for the sample of lEttori et alj j2010h . 
The best fitting slope is a = 0.48 + 0.09, substantially larger than 
the the oretical expectation and of the original finding of Etto ri et al.l 
(2010). This difference is likely due to the (1 + z) factor correction 
that they perform and to the different fitting techniques adopted. It 
should be noted that, contrary to the other works I consider here, 
the mass and concentrations of lEttori et alj {2010) refer to the dark 
matter only distribution rather than to the total mass. It is something 
to keep in mind for the interpretation of the results, however it does 
not change my conclusions, as discussed in Section[5] 



2.1.4 \Osuriet al\\20lA) 

Very recently, these authors used combined weak and strong lens- 
ing information in order to estimate the mass profiles of a sample 
of 28 galaxy clusters selected for having bright gravitational arcs 
in the SDSS cluster sample. The resulting catalog covers a rather 
large range in virial mass, 3 x 10 13 /r'M < m < 10 15 /r'M , 
and the redshift range where gravitational lensing is typically more 
efficie nt for sources at z s ~ 2, that is 0.3 < z < 0.6. 

lOguri et al] d201ll) acknowledge the recently established lack 
of redshift evolution for the concentration a t a given mass (at 
least for the lar gest objects, IGao et al.1 120081 ; IPrada et al.1 l201ll : 
IZhao et alj [2009). and therefore do not correct the measured con- 
centrations according to any redshift dependence. They find a sur- 
prisingly high slope for the concentration-mass relation, a ~ 0.6, 
substantially higher than that inferred from numerical simulations. 
Very interestingly, the authors propose a modification of the the- 
oretical concentration-mass relation that takes into account the 
strong lensing bias, namely the fact that strong lensing systems 
have on average a larger intrinsic concentration, and are preferen- 
tially elongated along the line of sight with respect to the general 
cluster population. Since this bias is more marked for low-mass 
systems, this increases the slope of the concentration-mass relation 
up to a ~ 0.2, sti l l muc h lower than found observationally. 

lOguri et al] d20 1 lh identify the problem as lying in the low- 
mass systems, as the high mass ones seem to be quite compatibl e 
with their model. This is a similar conclusion to lBuote et aljd 2007 ). 



In the middle right panel of FigureQ]l show the lOguri et al.l J201 lb 
sample, compared with the theoretical expectations. The observed 
slope that I find is a = 0.67 ± 0.13, in good consistency with the 
original estimate of the authors, and the highest slope of all the 
samples considered here. This latter fact can be ascribed to this 
being the only purely strong lensing selected sample, in which the 
lensing bias must play a role (although not sufficient to reconcile 
the observed slope with theoretical predictions). 

2.7.5 \Schmidt & AlieA\200% 

In this work the authors used Chandra X-ray observations of 34 
massive and relaxed galaxy clusters in order to infer the respec- 
tive mass distributions. The mass range they explored is relatively 
narrow, being limited to rich groups and massive clusters, m > 
3 x 10 14 /r'M . The studied systems display a variety of redshifts, 
ranging from almost local objects up to very distant ones (z ~ 0.7). 
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Figure 1. The observed relation between the concentration and the mass for groups and clusters of galaxies. Each panel refers to a different catalog of objects, 
compiled by the labeled authors. The blue long-dashed line represents the concentration-mass relation predicted by Gao et al (2008) at z = 0, while the solid 
magenta line is the best fit to the observed points, with labeled slope. The dotted magenta lines show the uncertainty on the best fitting slope. 
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ISchmidt & Allenl d2007l) managed to separate the dark matter halo 
mass profile from the distribution of the luminous matter compo- 
nents, however in order to compare with other studies, they also 
presented results of fitting the total matter distribution with a single 
NFW profile, which I adopted here. 

They compared their observational results with a power-law 
in mass having parameters estimated from numerical cosmological 
simulations, and in particular having a normalization that depends 
on redshift as oc (1 + zT l ■ Their results show that this model sub- 
stantially underestimates the concentration of low-mass systems, 
thus providing an overall poor fit to the observed concentration- 
mass relation. On the other hand, by relaxing the redshift depen- 
dence of the concentration normalization and allowing the slope of 
the concentration-mass relation to be free as well, they found a sig- 
nificantly better fit, having both a slope substantially higher than 
predicted by theoretical models (a = 0.45), and a remarkably small 
redshift dependence oc (1 + z )°- 71±052 (in terestingly in agreement 
with the simulations o f iDuffv et alj|2008h . All in all, this confirms 
the steep slope of the concentration-mass relation found in other 
works, and the mild redshift dependence of the concentration at 
fixed mass derived from the most recent numeri cal studies. In the 
bottom left panel of FigureQ]l show the results of lSchmidt & Allenl 
d2007t) . together with their best fitting power law. I obtain a slope 
of oc — 0.36 ± 0.14, consistent with the original estimate of the au- 
thors. This value is much larger than the expectations of theoretical 
investigations. 

2.1.6 \Woitak & Loka3H20ldl) 

In this last work, the authors employed kinematic data for a sam- 
ple of 41 nearby (z < 0.1) and relaxed galaxy clusters with 
m > 10 14 /j~'M , in order to constrain their mass profiles. The 
relaxation of a system was evaluated on the basis of X-ray ob- 
servations and the line-of-sight velocity structure, specifically, all 
clusters having strongly asymmetric or substructured X-ray maps, 
as well as bimodal velocity distributions were discarded from the 
sample. 

Contrary to the pr evious works summarized in this Section, 
Woj ta k & Lokasl feOld) found a slope of the concentration-mass re- 
lation consistent with flatness. On the other hand, the measured nor- 
malization is somewhat higher than predicted by theoretical mod- 
els, resulting in a marginal inconsistency with a WMAP-based stan- 
dard cosmo logy. In the bot t om ri ght panel of Figure Q] I report the 
results of Woit ak & Lokasl d2010b . I also found consistency with 
a flat concentration-mass relation (a = 0.03 ± 0.12), however the 
observed data do not seem to be in strong disagreement with the 
iGao et all d2008h prescription, although the measured concentra- 
tions seem to be somewhat systematically larger than theory would 
suggest. 



2.2 Implications 

The scrutiny of previous works presented in Subsection 12. 1 1 above 
shows that in general the concentration-mass relation of observed 
clusters and groups of galaxies is substantially steeper than any the- 
oretical prediction. In the two cases in which the measu red slope is 
compatible with the predictio n of the | Gao et al] d200 8) prescrip- 
tion, namely for t he sam ples of lComerford & Natarajan (200^) and 
IWoitak& Lokasl d201 Oh . the normalization results too high to be in 
agreement with a ACDM cosmological model having <x 8 ~ 0.8, as 
suggested by the latest WMAP data. Interestingly, this conclusion 



remains true quite irrespectively of the cluster selection function 
and analysis procedure, holding for either strong lensing clusters 
and for X-ray selected ones. In the former case one actually ex- 
pects some level of discrepancy with respect to theory, since strong 
lensing clusters are usually biased to have a larger inferred concen- 
tration than the average of the entire populati on, and the bias tends 
to be stronger for low er mass systems (see iHennawi et afl l2007l : 
iMeneghetti et alj|20ich . In their recent work. lOguri et alj d2oTTr> 
proposed a simple scheme in order to account for this lensing bias, 
finding a somewhat larger normalization of the concentration-mass 
relation and a slope that however can be increased up to a ^ 0.20 
at the most. While this would pro bably be a better fit to the data 
of lComerford & Nataraianl d2007l) . it is way insufficient to match 
their own data. Plus, X-ray selected clusters and kinematic stud- 
ies are not subject to any lensing bias, and still often show similar 
discrepancies. 

It should be noted that the observational works considered in 
this paper are by no means the only ones on the subject. There are 
plenty of them that are either older or making use of more lim- 
ited cluster samples. Despite their conclusions being less strong, 
many of these works hint at discrepanci es similar to those high - 
lighted above. An example is the w ork of | Gastaldello et alj d2007l) . 
whose data have been recycled in 
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:2007). I also re- 
2008) for a paper 



finding agreement with the theoretical concentration-mass relation. 

The systematicity with which observed low-mass systems are 
over-concentrated with respect to the theoretical expectations calls 
for some mechanism that is ignored in ;i-body simulations and that 
would be more effective at group scale rather than cluster scale. 
In this work I investigated the most obvious ingredient, that is the 
physics of baryons. The presence of gas and stars in real struc- 
tures affects estimates of concentrations in a twofold way: i) ob- 
servations often infer the total mass distribution of galaxy systems, 
which includes dark matter, hot diffuse gas, and stars. Numerical 
simulations on which theoretical models are based instead follow 
only the evolution of dark matter; ») the very cooling of gas in 
the central regions of dark matter halos imply a contraction of the 
dark matter component itself, an effect that is dubbed baryonic con- 
traction (see lGnedin et alfcoi ll and references therein). Given that 
gas cooling and star formation are expected to be more efficient in 
low-mass objects, both these effects should be more enhanced in 
those systems. It is my aim to make use of a simplified spherical 
model that takes into account these two factors, in order to see if 
the impact of baryons can indeed account for the observed steep- 
ness/normalization of the concentration-mass relation. 

It is worth to be noted that in several of the catalogs considered 
here there are a few objects having significantly lower concentra- 
tions than the theoretical ex pectations. The same has been f ound in 
the weak lensing analysis of Mandelbaum, S eljak. & Hiratald2008h . 
It is quite obvious that the gas cooling and star formation have only 
the effect of increasing the concentration of a structure, hence they 
cannot account for those data. Strong energy feedback from AGN 
as well as observational biases might be at work for those struc- 
tures, which I shall discuss better in Section[5] 



3 A SIMPLE SPHERICAL MODEL 

Baryons are present in two phases inside their host dark matter ha- 
los. The first phase is a hot diffuse gas component, which is re- 
sponsible for the emission of thermal X-ray radiation. This com- 
ponent dominates the baryon budget in rich groups and clusters of 
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Figure 2. The temperature-mass relation. Red circles with errorbars rep re- 
sent the results of the stacking analysis performed by Dai et alj J2010h in 
five different richness bins, ranging from poor galaxy groups to relatively 
rich galaxy clusters. The blue dashed li ne is the sca ling law obtained by 
combinin g the observed cr-T relation o f lXue&Wij 12000) with the m-o~ 
relation of lEvrard et all (2008) at z = 0. 



galaxies. The second phase is made by the cold gas and stars that 
constitute galaxies, and that result from the cooling and condensa- 
tion of the hot phase in compact regions. This phase is dominant 
in small groups and individual galaxies, since their shorter cooling 
time makes it easier to form massive stellar clumps. 



3.1 The baryon fraction 

The first fundamental ingredient for any self-consistent cosmic 
structure model lies in the fractions of hot gas and cold gas/stars 
that are present in their host dark matter halo, as a function of its 
mass and redshift. This issue is the focus of much active investiga- 
tion on the observational side, and a clear consensus has not been 
achieved yet. As of today it seems quite clear that the majority of 
baryons is outside bound structures, in the form of a diffuse and 
smoothly distributed gas component (Bregmarj2007h . However the 
actual fraction of baryons within dark matter halos and especially 
their partitioning between hot gas and stars is the subject of sub- 
stantial debate. In what follows, I shall introduce the prescription 
that has been adopted in this work, and then discuss at length how 
this compares with other wo rks in the s ame fi eld. 

In the recent paper bv lDaietal.1 J2010h the authors employ 
stacking techniques in order to derive the fraction of both hot gas 
and stars out to the virial radii of the host galaxy groups. With 
respect to other investigations on the same topic, this work has a 
twofold advantage from my point of view. First, it combines the 
stacked data for galaxy groups to pre-existing data on clusters as 
well as massive and dwarf galaxies, thus covering a very wide dy- 
namical range. The majority of other works focus instead on very 
limited mass ranges. Given that the concentration data I am con- 
sidering cover ~ 3 orders of mag nitude in mass, this is a desirable 
feature. Second. IDai et al. refer where possible their baryon 



fractions to the virial radii of structures (defined in the same way 
as done in the present work), while virtually all other works refer 
to radii encompassing overdensities of 500 or more with respect to 
the critical density of the Universe. This is another positive aspect, 
since many of the theoretical and observational works on concen- 
trations refer to virial quantities instead. Consequently, I adopted 
the fit to the baryon fractions measured in the work by iDai et al.1 

referring to that paper for additiona l details. 

There is one caveat though, namely IDai et alj l l201fj|) refer 
their results to the X-ray temperature or velocity dispersion of the 
structure at hand, thus forcing one to adopt scaling relations in or- 
der to convert them into a mass. As the authors themselves do, I 
convert the temperature into a velocit y disper sion by using the ex- 
perimentally calibrated relation of IXue & Wul 1 I2000I) , namely 



tr = 309 km s" 



kT 
1 keV 



(8) 



Then, by assuming that galaxies and dark matter particles are 
in dynamical equilibrium, I convert the velocity dispersion into 
a mass by a dopting the simulation-based scaling law derived by 
lEvrard et al.U2008l) . The combination of these two relations agrees 
reasonably well with the result of the stac ked analysis in the five 
richness bins adopted bv lDaietalJ (20 id) themselves, that cover 
the group-cluster mass scale. This is reported in Figure [2] As can 
be seen, the resulting observed temperature-mass relation tends to 
become somewhat steeper toward low-mass structures than it is at 
cluster scale. I tried to steepen a bit the association for low-mass 
groups, finding that the result is a substantial drop in the baryon 
fraction for m < 10 13 h~ l M Q , a regime that is basically of no inter- 
est here. 

The baryon fraction resulting from the analysis of IDai et"al] 

combined with the scaling relations mentioned above is de- 
picted in Figure [3] As can be seen, the stellar fraction increases 
monotonically with mass up to m ~ 10 13 tr [ M B , and then it starts 
a rapid decrease. That is also the same mass scale at which the gas 
fraction becomes as important as the stellar fraction, surpassing the 
latter for larger masses and being subdominant for smaller ones. 
The total baryon fraction always grows with mass, however it flat- 
tens at the mass scale of groups and clusters of galaxies, reaching 
the universal baryon fraction only for super-clusters. To my knowl- 
edge there are only two other works that push the measurement of 
the stellar mass frac tion (or baryon fract ion in general) out t o the 
virial radius, that are lConrov et alj (2007b and lAndreonl J20ich . For 
this reason, these findings ar e the only ones that can possibly be 
compared directly to the fit of IDai et all < |2010|) . 



The results o f lConrovet ail (2007) are represented in Figure|3] 
by the yellow triangles with errorbars, while the results of lAndreorJ 
(2010) are shown as magenta stars. As can be seen, the latter data 
are in rather good agreement with the prescription adopted here, al- 
though the stellar mass fraction seems to evolve som ewhat more 
steeply with mass. The data of IConrov et al] d2007b also are in 
agreement with the prescription of Dai et alj < |2010|) . except for the 
highest mass poi nt, which is lower by a factor of ~ 4. I note that 
the data points bv lConrov et al.1 (2007) do not take into account the 
stellar mass contribution of satellite galaxies, hence it is likely that 
the total stellar mass fraction displayed by the rightmost yellow 
triangle be significantly underestimated. It should also be kept in 
mind that the different works mentioned above (and below) adopt 
different conversions between luminosity and stellar mass, which 
can easily lead to ~ 50 - 100% discrepancies in the stellar frac- 
tion. Since I am only interested to gauge the possible impact of 
baryons on the concentration-mass relation, performing a highly 
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Figure 3. Th e baryon fraction as a function of total mass for z = 0, as given 
by the fit of iDai et alj feOld) . Each line style and color refer to a differ- 
ent matter component, as labeled in the plot. The horizontal cyan dotted 
line represents the universal baryon fraction given by the WMAP— 7 cos- 
mological parameters adopted in this work. The yellow t riangle s represent 
the stellar fractions measured at z = by Conroy et al. 1 2007 1 ) , whil e the 
magenta stars show the same quantity as measured bv lAndreonl 12010I) . Er- 
rorbars are not shown here since they are of the same size of the points. The 
thin lines show the gas and stellar mass fractions obtained by changing one 
of the parameters of the model, as explained in detail in the text. 



detailed analysis of these co nversions is bey ond the scope of the 
paper. I will hence adopt the IDai et al] d2010l) prescription as the 
fiducial one here, recalling that significant fluctuations around it 
are allowed. 

Several additional works studied the issue of the baryonic con- 
tent within cosmic structures, both from the points of view of the 
stars and of the hot diffuse gas. As mentioned above, all these works 
refer to substantially smaller radii than the virial ra dius, hence in or - 
der to make a fair comparison with the model of IDai etHll J2010h 
one should convert the masses (assuming specific profiles) and ex- 
trapolate somehow the baryon fractions. This latter operation is 
highly uncertain, as the distribution of baryons in the outskirts of 
cluste rs (dumpiness, depletion, etc., see for instance [Eckert etldl 
1201 lh is still very poorly understood. Therefore I do not attempt 
here such a procedure, rather I discuss only the most generic fea- 
tures of the baryon fraction, such as the slope as a function of mass, 
the scale of turn-around for the stellar mass fraction, etc. 

As already alluded to, the results on the stellar mass fraction of 
lAndreorj fcoiO) are in fair agreement with the model adopted here. 
The author also combined these findin gs w ith gas mass fra ctions 
measurements bv lVikhlinin et alj d2006j) and lSun et alj d200l^) . and 
upon correction for the different radial coverage adopted they con- 
cluded that the stellar and gas mass fractions are equal at m ~ 
10 13 /r'M Q , in agreement with Figure [3] Also, the total baryon 
fraction is a very shallow function of mass at cluster scales, again 
in agreement with the model adopted h ere. In a less recent X-ray- 
based analysis, Gast aldello et alj d2007l) found a significant decre- 
ment of the fraction of diffuse gas in galaxy groups as compared 



to cluster-scale objects, as can be seen in the fit of IDai et all d2010l) 
depicte d in Figure[3] 

In lGiodini et alj d2009h the authors studied the baryon fraction 
of a large sample of X-ray selected groups. They concluded that, 
when extending their mass range by considering previously ana- 
lyzed galaxy clusters, the stellar fraction scales with the mass with 

a power-law exponent of -0.37, compatible with the 0.4 that I 

find at the high-mass end for the model adopted here. In the same 
work, the gas fraction is also studied, based on previously existing 
compilations. The result is that the gas fraction increases with in- 
creasing mass , although the slope is quite modest (~ 0.2). Finally, 
iGiodini et alj d2009h also found a significant increment of the total 
baryon fraction as a function of mass at cluster scales, in agree- 
ment with Figure [3] of this paper. However, their slope is only of 
0.09, and they admittedly do not include the contribution from the 
Intra-Cluster Light (ICL), which mig ht flatten it further. The latter 
contrib ution has been studied by iGonzalez. Zaritskv, & Zabludoflj 
(2007). These authors found that the stellar mass fra ction changes 
somew hat more markedly w ith mass than found by IGiodini et al.1 
(2009) and lDaiet all (2010). while the change in diffuse gas frac- 
tion at cluster mass scales is compatible. Interestingly, the slope 
for the change of the gas mass fr action with mass (~ 0.21) is also 
consis tent with the one found by Pratt et al. (2009). Most remark- 
ably, IGonzalez. Zaritskv. & Zabludofij ( I2007D found that the total 
baryon fraction is constant for massive groups/clusters, advocating 
the impact of the ICL for t his discrepancy w ith other works. 

In a different study, Balogh et al. (2007) investigated the stel- 
lar mass fraction for a sample of galaxy groups, extending their 
analysis out to the virial radius. They found a definite decrement 
in the latter as a function of halo mass, however the uncertainties 
displayed in their results are very large (ranging from a factor of 5 
up to one order of magnitude), and for this reason I did not consider 
them in Figure[3] Nevertheless, those results seem to point toward 
a steeper decrease of the stellar f raction as compared to the model 
adopte d here, consistently with IGonzalez. Zaritskv. & Zabludoflj 
d2007l) . 

I conclude this Subsection by considering two very recent 
works that return results very different from the bulk of previous 
investigations w i th resp ect to the stellar mass fraction. Specifically, 
iLeauthaud et al.1 d201ll) found a content of stars for a sample of X- 
ray galaxy groups substantially lower than previous estimates (fac- 
tors betw een 2 and 5). Their results are m arginally consistent with 
those of iBehroozi. Conrov~ & Wechslei] d2010l) . that return even 
lower values for the stellar mass f raction. Interestingly the highest 
mass point of IConrov et alj d2007t) (which misses the contribution 
of satellites) is compatible with these works. Anot her difference be- 
tween these results and those derived by, e.g., the lDai et al.1 d20ld) 
analysis lies in the fact that the turn-around of the stellar mass frac- 
tion is at substantially lower masses, ~ 5 x 10" /r'M . It is not yet 
clear what is the main reaso n behind these largely discrepant re- 
sults. lLeauthaud et alJd201 lT) advocate inaccurate estimators for the 
stellar masses in previous works, however some problem might as 
well reside in the abundance matching technique adopted by these 
authors. 

In summary, excluding the works by I Leauthaud et al.1 d201lh 
and IBehroozi. Conroy. & Wechsletj d2010l) . the baryon fraction 
derived by the data of IDai et al.1 d20ld) in combination with 
suitable scaling relations is in good qualitative agreement with 
other works on the topic. Possible indic ations of a steeper 
slope for the stellar content (as foun d by Balogh et al. 20071 
IGonzalez. Zaritskv. & Zabludoflj 120071 and lAndreonl l201fj|) are 
only marginal and likely irrelevant for our conclusions. A possi- 
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Figure 4. The density profiles of different matter components in a system with total mass m = 10 13 /i~'M Q (left panel) and one with total mass m = 10 15 h~ l M Q 
(right panel). Each line refers to a different component, as labeled, while the green triangles show the total mass density profiles. 



ble di screpancy is the finding of iGonzalez. Zaritskv. & Zab ludoff 
J2007I) that the total baryon fraction should be constant wi th mass at 
cluster scales. Although the slight increase found by, e.g.. lAndreonl 
is consistent with those indications, this has not been 
stressed by other authors, and it might also be due to the different 
radii at which mass fractions are referred. In order to accommo- 
date for a substantial reduction in t he amount of star s in hig h-mass 
structures, as advocated recently by Leaut haud et alj J201 1 ). I mul- 
tiplied by 2 the slope of the correlation between the stella r to ga s 
mass ratio and the virial temperature presented in lDai et ah 
The resulting gas and star mass fractions ar e shown as thin lines in 
Figure [5] While keeping the original fit o f iDai et alJ d20ld) as the 
fiducial model, I will later show results for this latter configuration 
as well. 



3.2 Matter distribution within spherical systems 

Given the mass fractions of the different matter components that 
are present inside any given object, it remains to be established 
how these components are effectively distributed, i.e., their density 
profiles. Before looking into that, let me set up a little of notation. I 
shall indicate with f g (m, z) the fraction of hot diffuse gas present in 
a structure with total mass m placed at redshift z- Likewise, /*(m, z) 
and f.(m, z) represent the fractions in stars and dark matter for the 
same structure. It is obvious that /. (m, z) + f g (m, z) + f*(m, z) = 1 . 

Now back to the density profiles, the dark matter distribution 
prior to the contraction due to baryonic cooling shall be represented 
by a NFW d ensity profile, w ith the concentration-mass relation 
suggested by iGao et a I] d2008t) . The only caution to be exerted is 
that the mass of the dark matter halo (to which the concentration is 
attached) is not the total mass m of the structure at hand, rather it is 
f.(m,z) m. Also, it should be kept in mind that the cooling of gas 
and the formation of stars cause the dark matter mass distribution to 
become more compact than predicted by pure rc-body simulations, 
so that the final dark matter density profiles will not be NFWs any- 
more. I shall go back to this issue later below in Subsection l3.3l 



For the hot diffuse gas com ponent I adopted a y3-profile 
dCavaliere & Fusc o-Femiano 1976), which reads 



(1 +y2)W2 ' 



(9) 



where y = r/r c , This cored profile has formally three parameters, 
the outer slope /?, the core radius r c , and the corresponding core 
density p c . The latter however acts just as a normalization, whose 
value is determined by matching the total gas mass within the struc- 
ture with the gas fraction defined in the previous Subsection [TT] I 
shall assume as a characteristic object size the virial radius of the 
dark matter component prior to baryonic contraction, labeled as R. . 
This is very similar to the virial radius of the structure as a whole, 
having the advantage of being independent on the distribution of 
baryons. Hence, the following relation for the hot gas is enforced, 



4 , /3 3 5 

3*rPc (r c y.) 2F1 1 -, -p, -; -yt } = f g (m, z) m 



(10) 



In the previous equation y, = R./r c , while 2 Fi(a,b,c;x) is 
the Gauss hypergeometric function. This leaves the two free pa- 
rameters r c and p only. For the latter I adopted the reference 
value p = 2/3, which is in good agreem ent with the average 
outer gas slopes in observed X-ray cluster s ( Neumann & Arnaud 
Il999l : lOta & Mitsudal [20041: IVikhlinin et alj|2006t ICroston etaLl 
2008). Some authors dCroston et al.l I2OO8) find a dependence of 
this outer slope on the gas temperature, in that hotter clusters 
would have a steeper slope as compared to cooler ones. How- 
ever I did not include such a dependence, since the actual de- 
tails of the gas distribution are likely to play a marginal role in 
the inner part of struct ures, which is of interest her e. For the 
core radius r c I refer to[Ma kino. Sas aki. & Sutol d 19981) (see also 
ICapelo. Nataraian. & Coppl 2010h . where the authors have shown 
that a isothermal gas in hydrostatic equilibrium within a NFW dark 
matter halo is well represented by a /3-model with a core radius of 
about one fifth of the scale radius. For a typical halo concentration 
this corresponds roughly to r c = 0.05 R„ which I adopted through- 
out. Please note that for my choice of p = 2/3, Eq. ( |10t simplifies 
to 
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Figure 5. The contracted dark matter profiles divided by the original NFW profile, for a structure of total mass m = 10 13 h~ 1 M e (left panel) and m = 
10 15 h~ l Ms (right panel). The blue short-dashed lin e represents the stan dard adiabatic contraction, while the red shaded region enclosed by the two long- 
dashed lines shows the modified contraction model of Gnedin et al. 1201 ll) . Note that the radial distance from the center of the halo is expressed in units of the 
virial radius of the original NFW halo. 



4np c rl [y. - arctg (y.)] = f g (m,z) m . 



(11) 



For the stellar c omponent I assumed it to be radially dis- 
tributed according to a |jaffeNl983l) profile, that reads 



P*(0 



Po 



(12) 



where f = r/r . This equals at assuming all the stellar contribu- 
tion to be concentrated in the Brightest Cluster Galaxy (BCG) of 
clusters and groups, ignoring the presence of other galaxies. This 
is a hypothesis that can certainly be lifted in future more detailed 
studies. For the moment I am interested only in keeping the model 
as simple as possible, and in providing a proof of concept regard- 
ing the impact of this stellar component on the concentration-mass 
relation. In Section [5] I elaborate more on this approximation and 
argue that a redistribution of stars at large radii is mostly degener- 
ate with the as s umed stellar mass fraction. 

The lJaffel dl983l) profile has an isothermal slope of -2 at small 
radii, in agreement with the profile observed for stellar do minated 
systems such as elliptical galaxies dKoopmans et alj|2009h . As be- 
fore, the normalization of the profile is set by the stellar mass frac- 
tion, that is 



4Wo 1 



1 



/* (m, z) m , 



(13) 



where = R./r . The remaining free parameter r , which equals 
the half-mass radius for the stellar distribution, should be a few per- 
cent of the virial radius of the dark matter halo. I adopt ro = 0.02 R. 
henceforth, however results obtained by adopting different values 
for this free parameter will be shown. For a massive galaxy cluster 
this choice translates into rp ^ 20 /T 1 kpc, which is a ty pical effec- 
tive radius dSand. Treu, & Ellisll2002l ; ISand et alj|200a) . Because I 
assume that all stellar mass is included in this central clump, the 
value of ro could possibly be larger than this by up to ~ 60%. I 
refer to Figure[8]below to show that the impact of such a change is 
modest. Note that, since r <K R., the left-hand side of Eq. < | 1 3 1 > is 
actually very close to its asymptotic value Anpor^. 



In Figure|4]l show the density profiles of the three matter com- 
ponents inside two systems of different masses, one typical of gi- 
ant ellipticals/poor groups and the other characteristic of massive 
galaxy clusters. The different relevance of the various components 
can be clearly seen, with stars becoming dominant only at a few 
percent of the virial radius, while the hot diffuse gas remaining al- 
ways subdominant with respect to the dark matter component. Also 
note how the gas component is much more important than the stel- 
lar component in the outskirts of structures, more so for more mas- 
sive objects. It is significant that the stellar contribution starts hav- 
ing an influence at radii corresponding to the smallest ones probed 
by X-ray or gravitational lensing studies, because this implies that 
the apparently small differences seen between the total density pro- 
files of structures with different masses are going to have a substan- 
tial impact. Turning the argument around, the choice of the smallest 
radius for profile fitting has a large relevance on the derived param- 
eters. 



3.3 Baryonic contraction 

I mentioned above that the cooling of gas and the formation of 
large amounts of stars in the very central regions of cosmic struc- 
tures is expected to drag along dark matter, thus effectively in- 
creasing the host halo concentration. There has been an intense 
debate in the literature about whether this contraction effect is 
real and to what extent it affects actual structures. While a vari- 
ety of observational works find evidence for baryon contraction 
dMinor & Kaplin ghal 2008; ISchulz et al l l20ld: ISonnenfeld et al.1 
1 20 111 ; see also additional references in Gnedin et al J 1 20 1 lh . sev- 
eral others recover cored dark matter profiles i Sand et all 1 20041 : 
iNewman et al]|2009l . l20T ]]), so that the observational picture is not 
entirely clear at the moment. In the present paper I assumed that 
baryon contraction does take place. This is also supported by a very 
recent comprehensive study showing the occurrence and amplitude 
of the bar yon contraction effe ct in a variety of simulations, pre- 
sented bv lGnedin etal] d201ll) . Although this study is not the last 
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word on the topic (see the discussion in Section [5] about the role 
of AGN feedback), it provides a convenient way to parametrize the 
impact of baryonic cooling as observed in numerical simulations. 

Let me start by describing th e standard adiaba t ic con traction 
model bv lBlumenthal et alj ( fl986l) ; lRvden & GunrJ dl987h . which 
assumes any structure as made of concentric spherical shells which 
contract isotropically by conserving angular momentum. Given 
this, the radius r\ that initially encloses a dark matter mass m.{r{) 
can be related to the final radius r 2 enclosing the same mass by 



[m.(ri) + m bi iOi)] rj = [m.Oi) + m ba (r 2 )] r 2 , 



(14) 



where mt,,i(/"i) is the total initial baryon mass within r\, while 
m b,ii r i) ls me total final baryon mass within r 2 . In order to deter- 
mine the amount of dark matter contraction due to baryonic cool- 
ing it is hence necessary to establish initial and final mass pro- 
files for baryons. I assumed that there are no stars at some initial 
time in the dark matter halos, so that all baryonic matter is com- 
posed by hot diffuse gas, rab,i(n) = ra g ,i(n). The gas is still dis- 
tributed as described in the previous Subsection 13.21 only its pro- 
file is differently normalized by replacing f g (m, z) in Eq. JlOt with 
f g (m,z) + fi,(m,z)- After the cooling and star formation has oc- 
curred, stars and gas are distributed in the structure as described 

above, m b , 2 (r 2 ) = m g , 2 (r 2 ) + m^C/A 

As it has been shown by iGnedin et alj ( l201ll) . the adia- 
batic contraction model always overestimates the actual amount 
of contraction measured in numerical simulations. Hence, the au- 
thors provided a simple modification to Eq. j 14b that more accu- 
rately capture the contraction due to baryonic cooling (see also 
IGnedin et alj2004h . This modification reads 



[m.(n) + m b ,i(ri)] n = [m.(n) + m bj2 (r 2 )] r 2 , 
where 



(15) 



(16) 



0.03 R. \0.03R. 

IGnedin et alj J20T l|) suggested adopting A = 1.6, while they could 
not find a unique value of w that would fit equally well the wide 
variety of simulations they were considering. This spread is proba- 
bly ascribable to the variety of formation histories of cosmic struc- 
tures, as well as to the different implementation of baryonic non- 
gravitational physics that different simulation works adopt. In what 
follows I shall conside r a range of values fo r w that covers the vari- 
ety of results found by IGnedin et alj J201 lh . that is 0.6 ^ w ^ 1.3. 
Please note that if A = 1 and w = 1 the modified contraction sce- 
nario folds back to the standard adiabatic one. 

In Figure [5] I show the effect of baryonic contraction on the 
dark matter profiles of structures with two different masses. As 
naively expected, the effect of contraction is milder for larger 
masses, due to the relatively smaller amount of baryonic cooling. 
Also, I observe that the parameter choice brac keting the actual con- 
traction measured in various simulations by IGnedin et al.1 d201ll) 
provides generically smaller dark matter accumulation in the cen- 
tral regions than the adiabatic contraction model, as indeed was 
found by the authors themselves. 

Although not visible in Figure [5] the mass contra ction actu- 
ally becomes an expansion for the IGnedin et al] j201lh parameter 
choice A = 1.6 and w = 1.3 at r ~ 0.3 R.. The reason for this 
is that at this radius f ~ R., and the mass profiles are normalized 
so that initial baryonic mass and the final baryonic mass are equal 
at the virial radius of the original dark matter halo. Since however 
the final baryonic mass decreases more steeply than the initial one, 
outside the virial radius the former becomes smaller than the latter, 



thus effectively leading to a slight expansion. Finally, it should be 
noticed that the differences between the standard adiabatic contrac- 
tion scenario and th e modified baryonic contraction advocated by 
IGnedin etai]j20 111) are more pronounced at very small radii, while 
they turn out to be relatively mild at r > 0.01 R., the radial range 
that shall be considered in what follows. 



4 RESULTS 

I now turn to the main results of this work, namely describing how 
the mass and concentration of a galaxy system change with respect 
to the input values when fitting the total density profile, including 
the effect of baryonic contraction on dark matter. It is quite intu- 
itive that the result of this procedure would depend on the radial 
range used in the fitting, with fits extending to smaller radii re- 
turning arguably larger concentrations, due to the stronger effects 
of star formation and baryonic contraction in the inner regions of 
structures. In order to address this issue, I computed the total den- 
sity profiles, including baryonic contraction, for two objects hav- 
ing masses m = 10 13 lr l M e and m = 10 15 h~ [ M at n = 16 radii 
logarithmically spaced between 10~ 3 R. and R.. Then, I fitted NFW 
profiles to these density profiles, and compared the resulting masses 
and concentrations to their input values. I repeated this procedure 
n - 1 times, each time removing the smallest radius from the fit- 
ting procedure. The results, as a function of the minimum radius 
adopted in the fitting are displayed in Figure|6] 

Let me focus first on the behavior of the mass. Introducing 
baryons and the effect of their cooling on a structure, while still fit- 
ting the total density profile with a NFW function causes the mass 
to be slightly underestimated. This effect is very mild however, be- 
ing at most at the level of ~ 20% when the entire radial range 
[10~ 3 , 1] R. is considered. Additionally, I note that there is basi- 
cally no difference between different prescriptions for the baryonic 
contraction. On the other hand, the behavior of the concentration is 
much stronger. Including very small radii in the fitting procedure 
results in substantially overestimated concentration values, up to 
almost an order of magnitude. As expected, this effect is largely re- 
duced if only relatively large radii r > 0.05 R. are considered. On 
a related note, it is interesting to point out that in the observational 
study of iMandelbaum rSeliak. & Hiratal ( 120081) . who found a fair 
agreement with theoretical concentrations, a relatively large inner- 
most radius has been adopted. I also observe that the concentration 
overestimation is larger for smaller masses and when the standard 
adiabatic contraction recipe is adopted, because in both cases the 
impact of baryonic cooling and star formation is larger than for, 
respectively , larger masse s and the more realistic contraction pre- 
scription of IGnedin et al.1 ( 1201 ll) . For r > 0.01 R., the value that 
will be used in the remainder of this work, the difference between 
different baryonic contraction recipes is relatively limited, being at 
most of a factor of ~ 2 - 3. 

The results shown in Figure [6] are interesting in their own re- 
spect. In cases in which radial profiles can be measured over a wide 
range of radii, for instance via gravitational lensing, these results 
tell one what is the minimum radius one should consider in the fit- 
ting procedure in order to obtain meaningful values of the concen- 
tration. In cases where only a limited and relatively central radial 
range is accessible to observations, as for instance in X-ray studies, 
these results tell one how much the resulting concentration esti- 
mates are biased. Figure[6]also makes clear that in order to evaluate 
the impact of baryonic cooling on the concentration-mass relation, 
it is necessary to first specify a minimum radius for the fitting pro- 
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Figure 6. The masses and concentrations resulting from fitting an NFW profile to the total matter density profile, including baryonic contraction. The trends 
as a function of the minimum radius used for the fitting are shown. In both panels the lower set of curves refers to the mass (if = m), while the upper set 
represents t he concentration ( f = c) in units of their true values, that is the total mass of the structure and the concentration of a NFW profile with that mass 
according to Gao et al. 12008). The blue short-dashed lines show the results of assu ming the standard ad iabatic contraction prescription, while the red shaded 
areas enclosed by the red long-dashed lines show the modified contraction model of Gn edin et alj 1201 ll) . Finally, the left panel refers to m truc = 10 13 h~ l M , 
while the right panel to mime = 10 15 /t'Mq. 
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Figure 7. The concentration-mass relation resulting by including the effect of baryonic cooling and star formation. The blue short-dashed line represents the 
result of assumin g the standar d adiabtic contraction model, while the red shaded area bracketed by the red lon g-dashed l i nes re fer to the modified contraction 
model of lGnedin et"ai] 1201 1 ). For reference, the blue long-dashed curve shows the original relation found bv lGao et alj 1200 8|) by using n-body simulations. 
The left panel refers to the fiducial model for the baryon fraction as a function of mass, while the right panel considers a substantially reduced star formation 
in high-mass structures (see Figure[5]and the related discussion in the text). 



cedure, as different choices will most definitely lead to different 
results. In what follows I shall focus on the choice r m j„ = 0.01 R.. 
For the most massive clusters this translate in r min ~ a few tens of 
h~ l kpc, comparable to, e.g., the innermost radius adopted in X-ray 
studies. This approach is obviously a simplification of the actual 
procedure, however I am interested here in deriving general con- 
clusions. In more detailed future studies it will be easily possible to 



adapt the selected radial range to the observation one is interested 
to match. 

I now turn to the main point of the present investigation, that 
is, understanding how the concentration-mass relation is affected 
by the presence and cooling of baryons. In Figure[7]l sh ow the re- 
lation between concentration and mass implied by the I Gao et all 
J2008I) prescription, as well as its modification resulting from the 
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impact of baryons, computed with different recipes for the con- 
traction of dark matter halos. Also, in the two panels of the same 
Figure I present results for the fiducial baryonic fraction and its 
modified version (see Subsection l3.lt that results in a substantially 
reduced stellar abundance in high-mass structure, both depicted in 
Figure [3] As expected from the previous analysis (see Figure [6} 
the concentration at a fixed mass is always increased by the pres- 
ence of baryons, more so for low mass objects. However, this trend 
is not monotonic, rather it reaches a broad maximum at around 
m ~ 10 13 /r'M and then decreases slowly for lower masses. The 
reason for this stands in the stellar mass fractions. As can be seen 
in Figure [3] the stellar fraction does reach a broad maximum at 
m ~ 10 13 h M Q as well. For masses smaller than that the rela- 
tive abundance of stars decreases, and hence so does the impact of 
baryonic cooling on the concentration of structures. 

Overall, the impact of baryonic cooling is of the expected 
form, resulting in both a steeper slope and higher normalization for 
the modified concentration-mass relation. However, when adopt- 
ing the fiducial stellar mass fraction, the effect on the slope is not 
as pronounced as it would be suggested by th e observed relation s 
summarized in Se ction [2l especially those of | Ettori et all i2010h . 
lOguri et al.l j201 lb . and lschmidt & AllerJ j2007b . WherTthis is re- 
placed with a stellar mass fraction that is steeper at the high mass 
end, as it seems to be suggested by certain studies (see the dis- 
cussion in Subsection I3.lt , the concentration-mass relation gets 
steeper as well. Although the new relation is not well represented 
by a power law, by a qualitative comparison with Figure [Tj it ap- 
pears evident that the model having substantially less stars at large 
masses is in better agreement with observations. I shall make a 
more quantitative assessment in this respect further below. 

In Figure[8]l show the impact of the scale radius r of the stel- 
lar density profile on the modified concentration-mass relation. As 
can be easily understood, an increment in this scale radius results in 
a less compact stellar distribution. As a consequence the total mass 
profile is less peaked and the estimated concentration for a given 
mass tends to be smaller. The effect is however relatively limited, 
in that doubling r results in a reduction of the best fit concentra- 
tion of only ~ 10%. One would expect the opposite effect when the 
stellar scale radius is decreased. However, while the overall mass 
profile does become more peaked, the effect of baryonic cooling 
also shifts at smaller radii, thus moving outside the radial range 
that is adopted for the NFW fitting. This latter effect counteracts 
the former, so that the resulting change in concentration is only 
very subtle. Given that changes due to the stellar distribution are 
comparable to or smaller than those induced by different prescrip- 
tions for the baryonic contraction of dark matter, in what follows I 
shall stick to the fiducial choice ro = 0.02 R., and only show what 
happen if the baryonic contraction recipe is changed. 

For completeness, in Figure [8] I also report two additional 
lines. The magenta one shows the result of fitting an NFW func- 
tion only to the dark matter distribution after contraction, rather 
than to the total mass distribution. The radial range is the same 
adopted above, the baryonic contraction is adiabatic, and the stellar 
mass fraction is the fiducial one. The resulting concentration-mass 
relation is rather similar in shape to the previous one, however it 
displays a substantially lower normalization. This is expected be- 
cause the contribution of stars, which makes the overall density 
profile substantially peaked, is now ignored. The increment in con- 
centration at a given mass is now driven only by the contraction of 
the dark matter density profile due to baryonic cooling. Analogous 
conclusions are reached if the modified stellar fraction is a dopted 
instead. I note that in the simulations of iDuffv et al.l |2010) bary- 
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Figure 8. The co ncentration-mass relation according to the prescription of 
lGaoetalJi2003) (blue long-dashed line) and according to the modified pre- 
scription presented here, adopting adiabatic contraction of dark matter ha- 
los. The blue short-dashed line refers to the fiducial value for the stellar 
scale radius, while the other lines show the results of modifications of the 
latter, as labeled in the plot. The magenta dot-long dashed line represents 
the relation obtained by fitting only the dark matter profile (after adiabatic 
contraction) with an NFW function, rather than the total mass profile. The 
cyan dotted line shows the result of fitting the total mass profile without any 
baryon contraction effect. 



onic physics causes the dark matter halo concentration to be over- 
estimated by ~ 10% at most on cluster scales, while in Figure [8] 
I find overestimates of ~ 20 - 40% on the same scales. This sug- 
gests either a baryonic contraction less efficient than suggested by 
the adiabatic model (as expected) or a stellar mass fraction signif- 
icantly lower than predicted by the fiducial models adopted here 
(see below for further evidence in the same direction). The cyan 
line in Figure [8] displays what happen when fitting an NFW func- 
tion to the total mass profile, but without any baryon contraction 
effect. Again, I adopted the same radial range used above and the 
fiducial stellar mass fraction. As can be seen, the inclusion of stars 
without any contraction of the dark matter profile and the inclusion 
of baryon contraction without any baryons have about the same ef- 
fect on inferred concentrations at all masses, although the former 
results in a somewhat shallower relation. Note that these two effects 
cannot be simply summed together. 

In Figure|9]l show the same observed concentration-mass rela- 
tions of Figure Q] compared with the modified theoretical prescrip- 
tion employed here. I show the results obtained with the reduced 
stellar fraction, as they provide a better agreement with observa- 
tions (see below) than the fiducial one. As can be seen in the Figure, 
the modified concentration model does provide a better qualitative 
fit to the o bserved relations. Exce ptions to this are given by the 
samples of iBuote et"aH b007b and lwoitak & LokasJ bold) . In the 
former case the modified prescription seems to somewhat overes- 
timate the observed concentrations at small masses, while in the 
latter case nothing conclusive can really be said, due to the large 
scatter in the data. 
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Figure 9. The observed relation between the concentration and the mass for groups and clusters of galaxies, as presented in Figure [T] The blue long-dashed 
line represents the relationship predicted by Gao et al (2008) at £ = 0. The other lines refer to the concentration-mass relation obtained by including the effect 
of baryonic cooling and star formation with different models for the contraction of dark matter (see the text and previous Figures for details). 
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In order to quantify if and at what level the new concentration- 
mass relation provides a better representation of the observed data, 
for each observational sample I computed the x~ quanti ty defined 
in Eg. (|7), by replacing c(m,z = | c ,a) first with the lGao et al.1 
(2008) prescription, and then with the modified recipe adopted 
here. I used the adiabatic contraction model for reference, keep- 
ing in mind that suitably tuning the baryonic contraction design 
and/or the typical size of the stellar distribution (see Figure[8} can 
lead to better results. As it turns out the goodn ess of fit g = Vy^ 
decreases in all cases, with the exception of the lEttori et alj J20101) 
cluster catalog. This means that in almost all cases the modified 
concentration-mass relation including the impact of baryonic cool- 
ing is a better representation of observations as compared to pure 
n-body prescriptions. T he improvement is only slight (g = 8.9 be- 
comes g = 8.6) for the Woita k & Lokasl feoich catalog, while in 
all other cases i t is more substant i al, the best improvement being 
obtained for the Schmidt & Allen ( 2007) catalog (from g = 6.7 to 
g = 3.9). For the lEttori et al.1 d2010h sample the goodness of fit is 
almost unchanged, going from g = 6.3 to g = 6.5, implying that the 
modified relation does not provide any better representation of that 
catalog. However, it should also be recalled that in this work the 
authors report the masses and concentrations of dark matter halos 
only, thus caution should be taken in the comparison. I get back to 
this issue in the next Section[5] 

If the modified concentration-mass relation used above is re- 
placed by its version resulting from the fiducial baryon fraction, 
hence with more stars in high-mass structures, the situation is re- 
versed. Mo re precisely, the goodnes s of fit worsens in all cases ex- 
cept for the lSchmidt & Allerj d2007l) sample, for which it goes from 
the original g = 6.7 down to g = 5.0. This confirms the qualitative 
conclusion reached above, according to which a substantially re- 
duced stellar fraction in high-mass structures constitutes a better fit 
to the data. 



5 DISCUSSION 

As mentioned sev eral times in the course of the paper, the work by 
lEttori et al] d2010h reports the concentrations and masses of dark 
matter halos rather than whole structures, hence a comparison with 
the modified concentrat ion-mass relation reporte d in Figure[9]is not 
entirely self-consistent. lsThrnidt & Allenl J2007h also reconstructed 
the dark matter density profiles separately, however in order to de- 
termine the concentration-mass relation they referred to the total 
mass distributions, hence making the previous discuss ion fair. As 
for the remaining pure X-ray study bv lBuote et al . ( 2007), th ere the 
authors employ results from p revious papers (e.g.. lHumphrev et al.1 
l200d l iGastaldello et al. 2007). However in these previous works the 
mass fitting was performed in several different ways, both including 
and excluding the stellar contribution, in the latter case effectively 
fitting the overall matter distribution. It is not c lear which measure - 
ments were actually adopted in the study of iBuote et alj d2007h . 
he nce for the sake of completeness I put it on the same footing as 
the lEttori et alj {2010) work for the purpose of this Section. 

By considering the modified concentration-mass relation re- 
sulting from fitting the dark matter distribution alone and again 
considering the modified stellar fractio n that prov i des a suppressed 
stellar abundance at high masses, the lEttori et ail d2010h goodness 
of fit improves sli ghtly (from g = 6.3 down to g = 5.7), while the 
IBuote et al. (2007) goodness of fit improves substantially, decreas- 
ing from g = 8.5 to g = 5.0. The reason for the latter improve- 
ment is clear. The modified concentration-mass relation presented 



in Figure[9]overestimates the observed data, while the original pre- 
scription based on «-body simulations underestimate them. On the 
other hand, by considering the dark matter component alone gives 
rise to a relation that is in between the two, and hence fits very 
well the data. Also in the lEttori et al. I d20ld) case, the modified 
concentration-mass relation obtained by considering the dark mat- 
ter distribution alone fits the observed data better than the original 
prescription and the modified one adopted in the previous Section. 
The improvement however is very modest: the observed relation- 
ship is simply too steep to obtain a good match with theoretical 
models. 

Although introducing the effect of baryonic cooling and con- 
densation does improve the agreement of theoretical models with 
observations, there is still substantial room for improvement. One 
detail in particular that seems difficult to reconcile are the con- 
ce ntration values at l arge m ass es, that at least in t he X-ray studies 
of ISchmidt & Allenl f2007h andlEttori et al.1 d2010h (but also in the 
lensing study bv lOguri et al JI2OIIT) seem to be significantly lower 
than expected even with respect to pure dark matter studies. On 
the modeling side, significant AGN activity could prevent the for- 
mation of a cuspy stellar profile, and at the same time would also 
somewhat flatten the dark matter distribution. The simulations of 
iDuffv et al.1 d2010h shown that this is indeed the case and, surpris- 
ingly, the effect seems to be slightly more pronounced for high- 
mass clusters rather than groups of galaxies, a trend that would go 
in the correct direction. 

At the same time, dark matter halos extracted from rc-body 
simulations also present a scatter in concentration at a fixed mass, 
which is around ~ 0.15 in logarithm or more, and which might 
bring the less extreme observed points into agreement with the- 
ory within the respective uncertainties. However the systematicty 
of the discrepancy might still be a concern. On the observational 
side it is well know that, even for the most relaxed clusters as those 
considered in the catalogs mentioned above, non-negligible sub- 
sonic gas bulk motions typically lead to an underestimate of the 
cluster mass based on the hydrostatic equilibrium hypothesis by 
~ 10 - 20% dAscasibar et al. [20031: iRasia. Tormen. & Moscardinil 
l2004lRasia et al.ll2006l . l2008l ; lMeneghetti et al.ll2010h . This could 
contribute to steepen the observed concentration-mass relation, 
however only if this effect is not active on low mass objects, an 
instance that is not straightfo rwardly verified. 

Overall in the studi es of Schmidt & Allenl d2007h . lEttori et al.1 
d2010h . and lOgurietai] d201lh . the observed concentration-mass 
relation seems to be even steeper than what suggested by the correc- 
tion due to star formation presented here. While for the latter work 
this additional discrepancy can be explained via a strong lensing 
bias, this does not apply for the X-ray-based studies. One possi- 
bility is that for these cluster samples the stellar mass fraction at 
masses ~ 10 14 /r'M Q is substantially larger than indicated by the 
simple model adopted here. However, this might be difficult to rec- 
oncile with the observations of the baryonic fraction summarized 
in Subsection l3.ll The above mentioned underestimation of X-ray 
masses based on hydrostatic equilibrium also goes in the correct di- 
rection, however it is difficult to see why this sho uld not be at work 
on the group scale. It should also be noted that lEttori et alj j201Ch 
excluded the central 50 kpc of each cluster in their X-ray analysis. 
Since this is a larger fraction of the virial radius for low-mass ob- 
jects than for high-mass ones, Figure[6]suggests that this additional 
selection effect would flatten the modified concentration-mass re- 
lation, thus going in the wrong direction. 

One point that deserves to be discussed further is the impact of 
uncertainties on the distribution of gas and stars within dark mat- 
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ter halos. It is well known th at gas profiles exhib it a wide range 
of densities in structure cores dArnaud et alJl201CI) . Core densities 
lower than expected for a isothermal /J-model would result in a re- 
duction of the overall structure concentration. Still, the contribu- 
tion of hot gas to the overall density profile is subdominant at all 
radii, and non-negligible only in the outer parts of the most mas- 
sive clusters (see Figure |4j- Therefore I do not expect this uncer- 
tainty to change any of the conclusions of this paper. For the stellar 
distribution the situation is different, since I assumed here that all 
stars are in a single clump at the center of the structure. Remov- 
ing ma ss from this stellar clump an d redistributing it as a NFW 
profile dLin. Mohr. & S tanford 2004) could be more realistic. The 
net effect would be a decrement in the estimated concentration, 
more pronounced for systems that are less dominated by a Bright- 
est Cluster Galaxy (BCG), that are massive clusters. This proce- 
dure and its effect are basically similar to the decrement in the stel- 
lar fraction at high masses described above, since satellite galaxies 
give a negligible mass contribution wit h respect to, and their distri- 
bution is shallower than dark matter dBudzvnski etai]|2012h . As a 
matter of fact we verified that a trend similar to the one depicted in 
the right panel of Figure [7] could also be obtained via a redistribu- 
tion of stars within the structure, rather than a change in the stellar 
mass fraction. 

Another issue that might be relevant for the present discussion 
is the redshift dependence of the concentration for a fixed mass, 
that I have been ignoring throughout the paper. It is relatively well 
established in numerical simulations that the concentration for the 
high mass clusters is rather constant in time, implying that such a 
dependence cannot be advocated in order to explain the undercon- 
centration that several studies indicate in that mass range. On the 
other hand such a redshift dependence, despite being substantially 
weaker than oc (1 + z) , could flatten some what the observed rela- 
tion when factored in. I verified this for the lEttori et al.l d2010h cat- 
alog, finding that this is not the case. The steepness of the observed 
relation is unchanged, and the goodness of fit does not change sig- 
nificantly. 

Finally, it is useful to spend a few words on the redshift evo- 
lution of the modified concentration-mass relation that has been 
used here. Going at redshifts larger than zero, the average profiles 
of the various matter components would tend to change. For the 
dark matter this change is relatively well defined, being driven only 
by the evolution of the concentration. For other matter components 
this is less well established, for instance there is an ongoing dis- 
cussion about the extent to which the stellar distributi on is more 
compact at high-z dTruiillo. Ferreras, & de La Rosal201 11) . Another 
ingredient that might change with cosmic time is the baryon frac- 
tion, although no convincing evidence has been found for a sig- 
nifican t evolution of bo t h the gas and stellar mass fractions f rom 
z ~ 1 lAllen et alfeobl ; lEttori et1ii1l2009l : lG7odini et alj|2009b . In 
any case, once one's favorite recipe for the redshift evolution of 
these ingredients is set, it is straightforward to obtain a modified 
concentration-mass relation at any redshift. 



6 SUMMARY AND CONCLUSIONS 

In the present work I investigated the concentration-mass relation 
inferred from six different observed cluster catalogs, and tried to 
produce a modified theoretical relation that would take into account 
the role of baryonic cooling and star formation. The main results 
can be summarized as follows. 

• In almost all cases the slope of the observed concentration- 



mass relation is substantially larger than theoretical predict ions. 
The exceptions are give n by fhe lComerford & Nataraianl J2007b and 
Woita k~& Lokasl d2010b catalogs. In the former the slope is compat- 
ible with theoretical expectations, but the normalization is signifi- 
cantly higher than for a WMAP-based standard cosmology. In the 
latter the concentration at a fixed mass displays a very large scatter. 

• The systematicty with which the slope/normalization of the 
observed concentration-mass relation exceeds the theoretical ex- 
pectation indicates that low-mass clusters and groups of galaxies 
are over-concentrated, while massive systems have the expected 
concentration, or are possibly underconcentrated. Baryonic cool- 
ing affects more low-mass systems, hence it is expected to be a 
possible explanation for this discrepancy. 

• A simple model in which every structure is spherical and com- 
posed by stars, diffuse gas, and dark matter undergoing a contrac- 
tion due to baryonic cooling does indeed result in a larger estimated 
concentration for a fixed mass. As expected the effect is more en- 
hanced at the low-mass end, with concentration values being in- 
creased up to a factor of 2 - 3 for light galaxy systems. 

• This simple model also straightforwardly permits to determine 
the bias introduced in the measurements of the mass and concen- 
tration of an object as a function of the radial range covered. This is 
going to be extremely useful for the correct interpretation of future 
observational results in light of theoretical expectations. 

• The overall behavior of the resulting new concentration-mass 
relation depends significantly on the details of the stellar mass frac- 
tion. Low stellar abundances (or, alternatively, more extended stel- 
lar distributions) in large systems seem to be favored because oth- 
erwise the concentrations of these systems would be substantially 
overestimated with respect to the observed values. This result can 
provide a contribution to sort out the ongoing controversy about the 
stellar mass fraction in galaxy clusters. 

• The modified concentration-mass relation provides a better fit 
to the observed relation than theoretically motivated prescriptions 
in all cases (including those in which a steepening of the relation 
was n ot deemed necessar y), with the exception of the sample stud- 
ied bv lEttori et ail J201dh . 

Although very simplified in its nature, the modified 
concentration-mass relation that was proposed in this work pro- 
vides a definite improvement over the standard prescriptions based 
on fits to n-body simulations. Since it straightforwardly allows to 
include one's favorite baryon fraction recipe or density profiles, 
it will prove extremely useful in the interpretation of future sta- 
tistical studies on the matter distribution in galaxy groups and 
clusters. Possible advances include the introduction of a distri- 
bution of profiles for both stars and dark matter, as well as the 
role of ellipticity and substructures (satellite galaxies) on the fit- 
ting of a spherically averaged shape. As another example of ad- 
ditional effects that have not been included here, se veral authors 
dRudd, Zentner. & Kravtsovll2008l ; IZemp et al.ll201 lb have shown 
that the presence of baryons causes a global rearrangement of the 
total matter density profile resulting in a somewhat increased con- 
centration on top of the standard baryonic contraction. On a related 
note, substantial and complementary effort should also come from 
numerical simulations implementi ng baryonic physic s, of which 
one first example has been given bv lDuffvetaT]d2010h . 

At the same time, it is important to enlarge the observed 
sample of clusters and groups with measured density profiles and 
reduced systematic uncertainties. The current situation shows a 
very large sample-to-sample variance, as well as substantial scatter 
within each individual sample. The study of large and uniformly 
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selected cluster catalogs for which selection biases are well un- 
derstood will be a key step toward a better understanding of the 
formation of structures in the Universe. 
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